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NOTATION 
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stretching  stiffness 

cross-sectional  area  of  ring 

bending-stretching  coupling 

stiffnesses 

bending  stiffnesses 

Young’s  modulus 

shear  modulus 

total  in-plane  depth  of 

ring 

bending  moment 
circumferential  force 
concentrated  force 
internal  pressure 
shear  force 
mean  radius  of  ring 
shear  stiffness 
circumferential  displacement 
circumferential  displacement 
of  middle  surface 
radial  displacement 
radial  displacement  of 
middle  surface 


circumferential  position 
on  middle  surface 
outward  normal  position 
coordinate 

linear  and  quadratic 
coefficients  in  expansion 
of  V  in  powers  of  z 
shear  strain 
transverse  shear  strain 
at  z  =  0 
normal  strain 
circumferential  normal 
strain  at  z  =  0 
circumferential  position 
angle 

bending  curvature  change 
Poisson's  ratio 

normal  stress 
shear  stress 

cubic  coefficient  in  expansion 
of  V  in  powers  of  z 
i(  )/3y 


y 


1.  INTRODUCTION 


Ring-type  structures  are  an  important  class  of  structural  elements  used 
as  flywheel  rims,  shell  rims,  load-calibration  rings,  elements  of  transducers, 
and  as  circumferential  reinforcements  for  cylindrical  shells.  Furthermore, 
the  theory  of  rings  can  be  considered  to  be  a  degenerate  version  of  shell 
theory. 

The  theory  of  in-plane  vibration  of  relatively  thick  circular  rings, 
including  transverse  shear  deformation  and  rotatory  inertia,  dates  back  to 
the  work  of  Bresse  (1859),  which  predates  the  classical  thin-ring  theory  of 
Hoppe  (1871)  as  well  as  Timoshenko's  shear-deformable  beam  theory  (1921). 

A  discussion  of  some  of  the  large  body  of  ensuing  literature  was  provided  by 
Haines  (1974).  The  pioneering  shell  theory  due  to  Hildebrand  et  al .  (1949) 
and  the  work  of  Benzley  et  al .  (1973)  should  also  be  mentioned.  An  alter¬ 
native  is  to  use  the  two-  or  three-dimensional  theory  of  elasticity,  such  as 
used  by  Federhofer  (1935)  and  Endo  (1972),  respectively. 

It  has  long  been  known  that  the  Bresse-Timoshenko  theory  suffers  from  a 
serious  flaw:  it  is  based  on  the  assumption  of  a  uniform  distribution  of 
transverse  shear  strain  through  the  thickness,  whereas  the  actual  distribution 
is  that  of  a  distorted  parabola.  Numerous  ad  hoc  approaches,  both  static  and 
dynamic,  have  been  suggested  for  determining  the  so-called  shear  correction 
factor  to  account  for  the  differences  between  these  stress  and  strain  distri¬ 
butions;  see,  for  instance,  Mindlin  and  Deresiewicz  (1954)  and  Cowper  (1966). 

Only  recently  was  an  entirely  new  theory  proposed  to  allow  a  more  real¬ 
istic  shear  strain  distribution;  it  is  due  to  Levinson  (1981).  This  theory 
and  its  companion  plate  theory  (Levinson,  1980), as  well  as  analogous  earlier 
theories  by  Ambartsumyan  (1970),  Reissner  (1975),  and  Schmidt  (1977),  were 
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derived  for  macroscopically  homogeneous  material  only. 

It  is  well  known  that  transverse  shear  deformation  plays  an  even  more 
important  role  in  structures  composed  of  fiber-reinforced  composite  materials 
than  in  those  of  homogeneous  material;  see,  for  example,  Tarnopol ' ski i  et  al. 
(1965).  Dynamic  and  static  derivations  of  the  shear  correction  factor  for  a 
laminate  were  presented  by  Yang  et  al.  (1966)  and  by  Whitney  (1973)  and  Bert 
(1973).  The  latter  work  was  recently  extended  to  the  case  of  beams  laminated 
of  bimodular  materials  (different  elastic  properties  in  tension  and  compres¬ 
sion)  by  Bert  and  Gordaninejad  (1983). 

Levinson's  theory  of  plates  was  recently  extended  to  laminated  plates  by 
Murthy  (1981);  see  also  Bert  (1983).  However,  the  Levinson  theory  cannot  be 
extended  to  ring-type  members  and  still  satisfy  the  requirements  of  zero  shear 
strain  at  the  intrados  and  extrados  of  the  ring.  This  is  the  motivation  for 
the  present  theory,  in  which  these  zero-shear-strain  requirements  are  satisfied. 

Most  of  the  theories  of  shear-deformable  ring  mechanics  (including  Bresse, 
1859  and  Haines,  1974)  are  what  might  be  termed  first-approximation  theories, 
analogous  to  Love's  first-approximation  shell  theory  (see  Kraus,  1967),  in 
which  it  is  assumed  in  the  strain-displacement  relations  that  z/R  <<  1.  This 
is  analogous  to  Winkler-Bach  curved  beam  theory  (see  Seely  and  Smith,  1952) 
but  in  contrast  to  the  shear-deformable  ring  theory  due  to  Kirkhope  (1977). 

2.  HYPOTHESES 

In  the  present  work,  the  following  simplifying  assumptions  are  made; 

1.  As  in  other  ring  theories,  thickness  normal  strain  is  neglected, 
i .e. ,  ez  *  0. 

2.  The  shear  strain  distribution  is  nonlinear  in  such  a  way  that  the 
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shear  strain  vanishes  at  the  intrados  and  extrados  of  the  ring.  This  is  in 
contrast  to  all  other  ring  or  curved-beam  theories  (except  elasticity  theory) 
known  to  the  present  investigator,  and  allows  accurate  predictions  without 
the  use  of  a  shear  correction  factor. 

3.  The  theory  is  a  first-approximation  one,  i.e.,  z/R  <<  1.  This  means 
that  the  theory  should  be  most  applicable  to  thin  rings  made  of  highly  shear- 
deformable  material  (such  as  fiber-reinforced  composite  material). 

4.  Displacements  are  assumed  to  be  sufficiently  small  that  linear  strain 
displacement  relations  are  adequate. 

5.  The  ring  cross  section  may  be  homogeneous  or  arbitrarily  laminated  of 
linear  cylindrical ly  orthotropic  elastic  material  (see  Fig.  1). 


3.  THEORETICAL  DERIVATION 

The  following  displacement  field  is  used  as  the  point  of  departure 
(Fig.  1): 

V(\z)  =  v(.)  +  zc(:)  +  z  -  £  ( •'• )  +  z3:() 

(1) 

w(e,z)  =  w( o ) 


The  expression  for  circumferential  displacement  V  contains  the  same  terms  used 
in  Levinson's  beam  theory  (1981)  plus  middle-surface  displacement  v  to  account 
for  bending-stretching  coupling  and  the  quadratic  term,  which  is  now  necessary. 
The  expression  for  radial  deflection  W  is  consistent  with  hypothesis  1. 

It  is  shown  in  the  Appendix  that  it  is  not  possible  to  use  the  displace¬ 
ment  field  (1)  in  conjunction  with  the  exact  linear  strain-displacement  rela¬ 
tions  and  still  satisfy  hypothesis  2.  Thus,  consistent  with  hypotheses  3  and 
4,  the  following  first-approximation  strain-displacement  relations  are  used. 

e9  “  £y  “  <W/R>  +  V,y  ;  £z  *  W,z 
*  W  -  (V/R)  +  V  _ 


Yez  “  Yyz 


(2) 


6 


where  (  )  =  *(  )/3y,  etc.  Substituting  the  displacement  field  (1)  into 

»y 

(2),  one  obtains  the  strain  field 


£y  =  *-v,y  +  (w/R)^  +  Za,y  +  zZ\y  +  z’\y  ;  Ez  =  0 

V  ,  =  [w  -  (v/R)  fa]  +  [26  -  (a/R)]z  +  [3:  -  (2/R)]z:  -  (?/R)z? 
yz  »y 

The  material  stress-strain  relations  for  each  layer  in  the  laminate,  in 
view  of  hypothesis  5,  are 


y 


(4) 


In  view  of  the  third  of  (3)  and  the  second  of  (4),  enforcement  of  vanish¬ 
ing  shear  strains  at  the  extrados  and  intrados,  i.e.,  Tyz(y,  t  j  h)  =  0,  requires 
that 

[w  v  -  (v/R)  +  0  +  [3:  -  (,/R)](h/2):  =  0 
,y  (5) 

[26  -  (  :,/R)]  -  ( :/R)(  h/2)  =0 

p 

Consistent  with  hypothesis  3,  one  can  neglect  (1/8) (h/R)  compared  with 
unity.  Within  this  approximation,  it  can  be  shown  that  equations  (5)  are 
satisfied  if  one  sets 


6  =  (a/2R)  -  (l/6R)[w  -  (v/R)  +  i] 

,y  (6) 

C  =  -  (4/3h2)[w  -  (v/R)  +  a] 

>y 

Thus,  the  final  expressions  for  the  non-zero  strain  distributions  can  be 
written  as 

=  l\y  +  (w/R)]  ♦  Za>y  +  (z2/6R)[2a>y  -  w>yy  +  (V,y/R)] 

-  (4/3)(z3/h2)[Wjyy  -  (v>y/R)  +  a>y]  (7) 

vyz  -  [1  -  (1/3) (z/R)  -  4(z/h)2  +(l/6)(z/R)2  +  (4/3)(z3/Rh2)] 

•[w  „  -  (v/R)  +  a]  -  (a/2)(z/R)2 

9J 


The  nonlinearity  of  the  strain  distributions  in  (7)  is  readily  apparent. 
As  is  customary,  the  stress  resultants  and  stress  couple  are  defined  as 
(see  Fig.  1) 


(N,M) 


j  (l,z)oy  da 
a 


(8) 


From  the  theory  of  elasticity  in  plane  polar  coordinates  (R  +  z,')»  the 
circumferential  and  radial  equations  of  equilibrium  for  any  arbitrary  point 
within  the  ring  are,  in  the  absence  of  body  forces, 

t  _  +  (-.  .  +  2t  )/(R  +  z)  =  0 
T,z,,  +  (R  +  z)or,z  +  -r  '  =  0 


Multiplying  the  first  of  (9)  by  (R  +  z),  integrating  with  respect  to  z, 
using  relations  (6),  and  noting  that  t  vanishes  at  the  intrados  and  extrados, 
one  has 

N  v  +  (Q/R)  »  0  (10) 

Similarly,  multiplying  the  first  of  (9)  by  (R  +  z)z  and  integrating  with 
respect  to  z,  one  obtains 

M  -  Q  -  0  (11) 

«y 


Finally,  integrating  the  second  of  (9)  with  respect  to  z,  one  finds 


Q  „  -  (N/R)  +  p,  -  0 

yj  1 


(12) 


Substituting  the  strains  (7)  into  the  stress-strain  relations  (4)  and 
thence  into  the  generalized  force  definitions  (8),  one  obtains  the  ring  con 
stitutive  relations  as  follows 


’  N  \  [~A  B'lYv  y  +  (w/R)  1 

>J  ’  L8  D  -It  ’  J 


+  a) 


.y 


(13) 
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Q  =  $[w  -  ( v/R)  +  x]  -  S'-x 


04) 


where  the  ring  stiffnesses  are  defined  as 
(A,B)  s 


(1  ,z)E  da  ;  (C,F)  e  f  0,z)[(z2/6R)  +  (4/3h2)z3]E  da 


(B,D) 


(z,z^)[1  +  (z/2R)]E  da  ;  S‘  e  (1/2) 


(z/ R)  G  da 


05) 


S  E 


[1  -  (z/3R)  -  4(z/h)2  +  (1/6) (z/R)2  +  (4/3) (z3/Rh2)]G  da 


Equations  (13)  and  (14)  can  be  written  in  more  compact  form  as  follows: 


(:}•[ 


A  B' 
B  D 


Q  =  -S'., 

Substituting  the  ring  constitutive  relations  (13)  and  (14)  into  the 
equilibrium  equations  (10)-(12),  one  obtains  the  following  displacement 
equations  of  equilibrium  in  matrix-differential  operator  form: 

L1 

L, 

c 

L- 


Here, 


06) 

07) 


n 

L12 

L13 

r 

v/R 

r°] 

21 

L22 

L23 

< 

* 

0 

31 

L32 

L33 

l  W/RJ 

08) 


n  £ 

(AR  +  C)dy2  -  (S/R) 

,  l12 

13  s 

-  CRdy3  +  (A  +  S)dy 

,  L21  e 

22  5 

(D-F)dy2  -  S  +  S' 

HI 

CO 

CM 

t31  5  -  [(A  +  S)  +(C/R)]dy  ,  L32  H  [S  -  S*  -  (B'/R)  +  (C/R)]dy  , 

L33  =  (SR  +  C)dy2  -  (A/R)  ,  dy  =  d(  )/dy  (19) 

4.  APPLICATION  ANO  COMPARISON  WITH  EXPERIMENTAL  RESULTS 
As  an  application  of  the  new  theory,  a  diametrically  loaded  thick  circu¬ 
lar  ring  (see  Fig.  2)  is  considered.  Application  of  statics  to  the  free-body 
diagram  shown  in  Fig.  3  yields  the  following  expressions  for  the  generalized 
forces : 

N  =  -  P  cos  9  ,  Q  s  -  P  sin  9  ,  M  *  -  M1  + PR  cos  9  (20) 

These  expressions  satisfy  the  equilibrium  equations  of  the  present  theory, 
equations  (10)-(12). 

Since  the  problem  is  statically  indeterminate  to  the  first  degree,  an 
additional  equation  relation  to  deformation  is  required.  From  the  ring  con¬ 
stitutive  relation,  equation  (13),  one  has 


a  =  [AM  -  BN  +  (AF  -  BC)(Q  /S)]/[AD  -  BB'  -  (AF  -  BC)(S'/S)] 
>y  »y 

Integrating  symbolically,  one  obtains 


a(e) 


9 

a  R  d0  +  a(0) 
»y 

0 


(21) 


(22) 


Inserting  the  appropriate  quantities  into  equation  (22)  and  enforcing 
the  symmetric  conditions 

o(0)  =  a(ir/2)  *  0  (23) 

one  finds  the  following  expression  for  the  statically  indeterminate  bending 


moment: 


M]  *  (2/ir)[R  +  (B/A)  -  (1/ARS)(AF  -  BC)]P 


(24) 
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It  is  noted  that  in  the  case  of  a  homogeneous  beam  in  thin  ring  theory 
(S  -*•  “),  equation  (24)  reduces  to  the  classical  value,  M-j  =  (2/t)PR. 

The  author  was  not  able  to  locate  any  experimental  stress  or  deflection 
data  for  a  circumferentially-wound,  fiber-reinforced  composite-material  ring, 
which  would  be  the  most  stringent  test  for  tbe  new  theory.  Hbwever,  some 
photoelastically  obtained  stress  data  were  reported  by  Srinath  and  Acharya 
(1954)  for  a  very  thick  ring  of  homogeneous,  isotropic  material.  Specific 
data  for  their  ring  are  listed  in  Table  1. 

For  the  present  special  case  of  a  homogeneous,  isotropic  ring,  the  ring 
stiffness  parameters  can  be  expressed  as  follows: 

A/E  =  h  ,  B/E  =  0  ,  B'/E  =  h3/24R  ,  C/E  =  B'/3E  ,  D/E  =  h3/12  , 

3  2  (25) 

F/E  =  D/5E  ,  S/E  ■  (l/3)h/(l +v)  ,  S'/E  =  (l/48)(hJ/R  )/(l  +  v) 

It  was  necessary  to  normalize  the  stiffnesses  with  respect  to  Young's  modulus 
(E),  since  the  actual  value  of  the  Young’s  modulus  was  not  reported  by  Srinath 
and  Acharya.  The  resulting  expressions  for  the  circumferential  normal  stress 
distribution  are  given  below. 

For  the  actual  value  of  shear  stiffness  (S): 

oy  *  -  525  +  16,015  z  +  6,548  z2  +  13,195  z3  (26) 

Here,  z  and  oy  have  units  of  . inches  and  psi,  respectively. 

Neglecting  transverse  shear  flexibility  (S  -*■  ®): 

oy  =  -  517  +16,104  z  +  6,315  z2  -  10,477  z3  (27) 

It  is  noted  that  inclusion  of  transverse  shear  deformation  results  in  a 
more  nonlinear  stress  distribution.  A  comparison  of  these  predicted  results 
with  the  reported  experimental  values  is  presented  in  Table  2.  Certainly 


the  results  for  the  new  theory  (actual  S)  are  encouraging.  Also  presented 
is  the  prediction  of  Winkler-Bach  ring  theory,  for  which  the  bending  stress 
is  given  by 


~b  =  (M/he) (e  +  z)/(R  +  z) 


where  e  is  the  distance  from  the  centroid  of  the  cross  section  to  the 
neutral  surface  and  is  given  by 


e  =  R  -  h[ln 


1  +  (l/2)(h/R)-i~1 
1  -  (1/2) (h/R)J 


(28) 


(29) 


In  the  present  case,  at  the  horizontal  cross  section  (:=0),  the  bending 
moment  is  given  by 

M  =  P(R-e)  -  M,  =  [1  -  (2/-)]P(R-e)  (30) 


The  total  normal  stress  at  -  =  0  is  given  by 

c  •  -  (P/h)  +  ,-b  (31) 

In  the  numerical  example  under  consideration,  equation  (31)  becomes 

2  =  -  437  +  18,835  (0.0099356  +  z)/ (1 .1875  +  2)  (32) 

The  numerical  values  are  presented  in  the  next  to  the  last  column  in  Table  2. 

There  is  not  a  great  deal  of  difference  among  all  four  stress  distri¬ 
butions  listed  in  Table  2.  However,  it  is  emphasized  that  the  material  is 
isotropic  with  v  =  0.35,  thus,  the  ratio  of  shear  modulus  to  elastic  modulus 
(G/E)  is  0.3704.  In  contrast,  for  graphite/epoxy  (T300/5208),  this  ratio 
is  only  0.04368,  according  to  the  carefully  correlated  experimental  results 
due  to  Knight  (1982).  Thus,  it  was  decided  to  repeat  the  calculation  of 
the  present  theory  for  a  hypothetical  ring  constructed  of  hoop-wound 


graphite/epoxy  and  having  the  same  geometry  and  loading  as  previously  di 
cussed.  The  resulting  equation,  instead  of  equation  (26),  is 

oy  =  -  587  +  15,204  z  +  8,163  z2  +  119,000  z3 

This  calculation  was  based  upon  the  following  material  properties 
reported  by  Knight  (1982): 

Hoop-di rection  elastic  modulus  20.1  x  10^  psi 

Transverse  shear  modulus  0.378  x  10^  psi 

The  resulting  stress  distributions  are  plotted  in  Fig.  4,  where  it 
is  noted  that  the  low  shear  modulus  of  the  graphite/epoxy  causes  much 
higher  maximum  stresses  than  are  present  in  the  isotropic  ring. 


5.  DISCUSSION 


In  the  development  of  the  equilibrium  equations,  (10)-(12),  the  inte¬ 
gration  method  believed  to  have  been  originated  by  Yang  (1965)  was  used. 
However,  this  approach  leads  to  equilibrium  equations  which  are  only  approxi¬ 
mate  for  the  present  case,  although  they  are  exact  for  the  cases  of  the  clas¬ 
sical  thin  and  shear  deformable  theories.  More  accurate,  yet  much  more  com¬ 
plicated,  equilibrium  equations  can  be  obtained  by  application  of  the  princi¬ 
ple  of  virtual  work,  such  as  used  by  Kraus  (1967)  for  thin-shell  theory.  How¬ 
ever,  in  the  present  case,  this  approach  gives  rise  to  five  more  generalized 
forces  than  the  theory  presented  here.  They  are  of  the  form 


,h/2 

i  (z'V  2  V  2Tyz’  z  Tyz ’ 
-h/2 


These  are  analogous  to  the  higher-order  generalized  forces  obtained  by  Tiffen 
and  Lowe  (1963)  and  Lo  et  al.  (1977a, b)  in  their  higher-order  plate  theories. 
See  also  the  treatise  by  Librescu  (1975).  The  present  approximation  is  more 
in  keeping  with  Levinson's  plate  and  beam  theories  (1980,  1981).  Thus,  it 
may  be  more  appropriate  to  classify  the  present  theory  as  a  more  accurate 
simplified  theory  rather  than  as  a  true  higher-order  theory.  It  should  be 
mentioned  that  Langhaar  (1965)  pointed  out  that  higher-approximation  theories 
are  not  necessarily  more  accurate.  He  cited  the  case  of  a  curved  cantilever 
beam  for  which  Love's  first-approximation  theory  (same  as  Winkler-Bach  curved- 
beam  theory)  yielded  the  statically  correct  force  quantities  while  his  second- 
approximation  theory  gave  results  which  violated  static  equilibrium  grossly. 
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6.  CONCLUSIONS 

The  new  ring  theory  presented  here  has  the  following  advantages: 

1.  In  contrast  to  other  shear-deformable  ring  theories  it  does  not 
require  separate  cid.  hoc  determination  of  a  shear  correction  factor. 

2.  In  contrast  to  other  fi rst-approximation  ring  theories,  it  allows  for 
a  nonlinear  distribution  of  circumferential  displacement  through  the  thickness. 

3.  Relative  to  two-  or  three-dimensional  elasticity  theory  and  to  higher- 
order  ring  theory,  it  is  much  simpler  and  much  more  amenable  to  engineering 
calculations. 
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APPENDIX 


The  exact  linear  strain-displacement  relations  in  shell  coordinates  are 

. 4  «  (W  ♦  V  )/(R  ♦  2)  ;  ,  *  W  , 

Z  ,2  { A- 1 ) 

y ..  *  V  ♦  (U  -  V)/(R  ♦  2) 

i  » 

Using  the  displacement  field  of  equations  (1)  in  equations  (A-l),  one 
obtains 

=  (R  ♦  2)*1  (w  ♦  V  *  2,  *  z2  *  23  ) 

*  *  *  • 

z  -  o  ;  .  z  =  (R  +  2)"1 (m  -  V  -  2  ,  -  2  -  2  )  ♦  .  ♦  2z.  ♦  32  :  (A-2) 

Enforcing  .  (  ,  ♦  ^  h)  *  0  thus  requires  that  the  following  four  equations 
be  satisfied: 

w  -  v  -  (h"/4;  =0 

( A  -  3 ) 

«  ♦  (h2/4)  *  0  ;  -  0  ;  .  ♦  (3h2/4)  =  0 

It  is  obvious  that  all  four  of  these  cannot  be  satisfied  except  in  the 
trivial  case,  «  a  .  *  *  0  and  w  =  v . 
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Table  1.  Conditions  for  experiments  of  Srinath 
and  Acharya  (19540 

_ Dimensional  Data _ Loading  and  Material  Data 

Ring  depth  h  =  0.375  in.  Loading  P  =  164  lb/in. 

Ring  mean  radius  R  =  1.1875  in. 

Ratio  h/R  =  0.3158  Poisson's  ratio  =  0.35* 

Estimated 


Table  2. 

Numerical  results  for 
normal  stress  (psi) 

distribution  of 

circumferential 

Present  Theory 

Winkler-Bach 

Theory 

z/h 

S  -» 00 

Actual  S 

Measured 

-0.500 

-3,246 

-3,371 

-3,781 

-3,550 

-0.233 

-1,871 

-1 ,882 

-1,765 

-1,720 

0 

-517 

-525 

-279 

Not  reported 

0.042 

?64 

-273 

-37 

-70 

0.291 

1,301 

1,322 

1,292 

1,300 

0.500 

2,656 

2,809 

2,268 

2,350 
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Fig.  1  .  Coordinates,  displacements,  forces  and 

moments  acting  on  a  differential  element 
of  a  ring. 


» 


POSITION  MEASURED  OUTWARD 
FROM  CENTROID,  IN. 


Fig.  4.  Normal-stress  distributions  for  isotropic  and 

graphite/epoxy  rings  loaded  in  diametral  compression 
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